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ABSIBACl 

Absolute rate constants for the reaction OH + H^CO have been measured 

over the temperature range 228-362K using the flash photolysis-resonance 

fluorescence technique. The results were independent of variations in [H^CO], 

total pressure [Ar] and flash intensity (i.e., initial [OH]). The rate 

constant was found to be invariant with temperature in thi3 range, the best 

-11 3 - 1-1 

representation being = (1.05 ± 0.11) x 10 cm molecule s where the 
error is two standard deviations. This result is compared with previous 
absolute and relative determinations of k^. The reaction is also discussed 
from a theoretical point of view. 


The reaction of the hydroxyl radical with formaldehyde is of interest in 

1 

high temperature combustion studies, in an explanation for CO formation in 

2 

the atmosphere of Jupiter, and in methane oxidation in the troposphere and 

stratosphere of the earth.- 5 In all these instances, H 2 C0 is a product of 

hydrocarbon oxidation. For the Jovian atmosphere, it is proposed that H p C0 is 

formed primarily in the reaction 0 + CH. . Loss processes include photolysis 

2a ' 5 

ana reaction with H and OH. . In the stratosphere, H 2 C0 is formed as a 

product of the low temperature oxidation of atmospheric methane. It is 

removed by photolysis and by reaction with OH, 0 and Ct. The reactions 0 + 

h„C0 and Cs. + H„C0 are relatively minor loss processes for formaldehyde at 
^ d 4 

most altitudes in the stratosphere. Rate constant measurements for these 

5 5 

reactions have been recently reported by us as well as other workers. ’ 

Reaction with OH is a major sink for H_C0 at all altitudes in the upper 
4 d 

atmosphere, but there have been no direct measurements of the absolute rate 
constant at the low temperatures prevailing there. 

The mechanism of the reaction has been usually assumed to be H atom 
abstraction 


OH + H 2 C0 + H 2 0 + HCO, (la) 

although there is little or no direct evidence for this. Recently, Horowitz 
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et al. have proposed the additional reaction channel 

OH + H 2 C0 * HCOOH + H. (1b) 

Additional work is required to test the validity of this suggestion. 

Room temperature measurements of the absolute rate constant for the 
reaction 


£ 





OH + H 2 C0 products ( 1 ) 

have been made using the discharge flow-mass spectrometric technique (DF-MS) 

8 9 

with OH in excess ' and by the flash photolysis-resonance fluorescence 

in d 

technique (FP-RF) with H 2 C0 in excess. Herron and Penzhorn 0 obtained a 

-11 ^ -1 -1 q 

lower limit k 1 > 0.7 x 10 cnr molecule s , Morris and Niki^ give k^ = 

(1.4 ± 0.3^) x 10” 11 cm 3 molecule” 1 s” 1 and Atkinson and Pitts 10 report k 1 = 

(0.94 ± 0.10) x 10” 11 cm 3 molecule” 1 s" 1 , all at 298 ± 2 K. There have been 

three relative rate constant determinations at room temperature. Morris and 

Niki, 11 using the DF-MS technique, measured OH + H„C0 relative to OH + C,H fi 

and obtained the result k, = (1.5 ± 0.1 C ) x 10 cm 3 molecule” s” 1 . Niki et 
12 1 5 

al. employed Fourier Transform Infra- red Spectroscopy ( FTIR) with OH + C ? H 2 . 

as the reference reaction and obtained k = ( 1 . 5 ± 0. 1 ) x 10 cnr molecule 
-113 1 

s . Smith measured OH + HgCO relative to OH + OH using the DF-MS technique 
and reports k^ = 0.65 x 10” 11 cm 3 molecule” 1 s” 1 . 

There have been only two variable temperature studies of the reaction, 
one absolute 10 and one relative. 13 Atkinson and Pitts 10 give k. = 1.25 x 

“1 1 3 _ -1 _ 1 ' 

10 exp(-88 ± 150/T)) car molecule s for the temperature interval 299 to 

13 

426 K while Smith's results from 268 to 334 K may be represented by k. = 6 x 

— 1 1 3 11 ' 

10 exp(-635 ± 250/T) cnr molecule s . This later study is the only one 

below room temperature but its usefulness in considerably reduced by the fact 

that the rate parameters of the reference reaction OH + OH are far more 

uncertain than those for OH + itself. 


An inspection of the available rate data thus shows disagreement by a 
factor of two at room temperature, a factor of five in the Arrhenius 
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pre-exponential factor and a factor of seven in the activation energy. All of 
these are well beyond any reasonable estimates of error. Further, as noted 
above, there are no reliable rate data at low temperatures. Because of these 
discrepancies and the need for low temperature rate data for atmospheric 
modeling, we have undertaken the present study. V»e report here absolute rate 
constant measurements over the temperature range 228-362 K. 

EXPERIMENTAL 

The present experiments were performed by the flash photolysis-resonance 

fluorescence (FP-RF) technique. The apparatus and specific techniques are 

14 i 

quite similar to those reported previously; ’ " however, we have replaced the 
reaction cell and vacuum housing assembly. 

The new cell is a black teflon coated brass cylinder whose internal 
dimensions are three inches in diameter and three inches in length. 

Fluoresence is observed coaxially in the center of one end wall of the 
cylinder. Honeycomo baffling occupies the other end of the cylinder. 

Resonance and flash lamp radiation enter the cell at right angles both to each 
other and to the detection port on the cylinder surface at about one and a 
quarter inches from the detection end wall. Optically flat windows or lenses 
(0.75 inch aperature) are compression sealed to the ceil through either Viton 
or silicone rubber "0" rings. The cell has entrance and exit tubes for 
premixed gas flows. Also entrance and exit tubes are attached for cooling 
gas. The cell has eight interconnected and equally spaced one sixteenth inch 
channels through the brass wall to provide a continuous path for the cooling 
gas. In addition it has eight equally spaced connected channels to 
accommmodate an insulated nicrome heating element. Temperature is measued by 
a chromel-alumel thermocouple which is positioned inside the cell near the 
fluorescence volume but out of the line of sight of the photomultiplier which 
detects the resonantly scattered photons. The cell assembly is placed inside 
a brass vacuum chamber which has optical access ports for detection, resonance 
light, and photolysis light admission. Thus, this newly constructed 
equipment, while retaining all of the functions of the previous design, has as 
its principal advantage significantly lower scattered background light 
compared to the earlier cell. 
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For the reported experiments, the source of OH radicals is flash 

photolysis of water. OH resonance radiation (* * 306 nm) is produced from a 

water saturated discharge in Ar at 17 torr. This radiation is highly 

collimated and slightly focused by a lens and collimator system before 

entering the cell. The fluorescent photons are collected and focused by a 

lens-collimator-interference filter (Ditric Optics, Inc., 310 nm, FWHM = 11 

nm) optical system onto an EMR 5^1 h photomultiplier tube. The fluorescent 

photons are recorded in repetitive flashes by a multichannel analyzer 

operating in the multiscaling mode. The photoflash light is also spectrally 

isolated by a Ditric Optics, Inc. interference filter centered at 157.5 nm 

with a skewed FWHM of -13 and +25 nm, the maximum transmittance being 28 %. 

This filtering technique also has the added advantage of decreasing the 

scattered light due to photo flash interference. The initial yield of OH from 

the photodecomposition of H 2 0 using this filter was 19% of that obtained using 

a LiF window (1 > 105 nm) . This factor reflects the spectrally integrated 

effects of wavelength dependent transmission, H 2 0 absorption, and photo flash 

intensity. We roughly estimate from a consideration of these factors and 

11 -3 

yields of H from H 2 0 photolysis, that initial [OH] is about 1 x 10 cm for 
these experiments. 

The experiments, and particularly the handling of H o C0, were carried out 

6 ^ 

in a similar way as described previously. We found no unusual behavior with 
HgO-HgCO-Ar mixtures as long as the flow rate of the premixed ternary mixture 
was sufficiently high. The present ceil is larger than that in our previous 
investigations, and we estimate that the contents of our new cell are 
completely replenished every 2-3 flashes. 

Experiments were carried out at four temperatures from 228-362 K with 
substantial variations of [h^CO] and total pressure. Due to the presence of 
the interference filter in the photoflash optical path and a practical upper 
limit on flash intensity, the variation of flash intensity is somewhat limited 
compared to that used in other studies from this laboratory. Over this 
limited range, however, no flash intensity effects were observed. Experiments 
were performed with and without added [H 2 C0] in order to assess the 
diffusional loss corrections for OH. The observed psuedo- first-order OH 
radical decay plots were strictly linear under both conditons, and the 
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diffusional corrections were typically less than 15$ except for some 
experiments at 362 K where this contribution was 35$. 


Argon (Matheson, 99*9995$) for mixture preparation and Argon (Goddard 
grade, water pumped) for the resonance lamp were both used without further 
purification. Deionized water was used as received in the discharge resonance 
lamp but was bulb to bulb distilled, the middle fraction being retained, for 
mixture preparation. Formaldehyde was prepared from paraformaldehyde as 
described previously. Mass spectral analysis showed the sample to be 99-7$ 
pure with water being the only measurable impurity. 

RESULTS 


Under the pseudo-first-order conditions employed in this study, with 
[H 2 C0] » [OH], the decay of OH radicals is given by 


S-n[OH] = - k 


observed 


t + «-n[OH] 

• o 


( 2 ) 


The observed pseudo-first-order decay constant is represented by 


observed 


= k. 


[H 2 C0] + k d 


(3) 


where k 1 is the bimolecular rate constant for reaction (1), and k d is the 
first-order rate constant for OH radical diffusion from the reaction volume 

viewed by the pnotomultiplier . k^served and k d are obtained by linear least 

squares methods from plots of the logarithm of accumulated counts against 
time, where accumulated counts are proportional to [OB], Typical results for 
the reaction at 228 K are shown in Figure 1. Values of k^ were calculated 
from k observed according to equation (3). The values obtained for are 
reported in Table 1 for the various conditions employed in this study. At 
each temperature, experiments were performed with variations in total 
pressure, [H 2 CO], and flash intensity. The errors are one standard deviation 
values as determined from repeat determinations at the indicated experimental 
condition. The flash energies reported are the effective flash energies 
entering the cell through the 157-5 nm interference filter. The results show 
that k,j is invariant with temperature over the range of 228 < T < 362 K. The 
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best representation for this temperature range is k„ = (1.05 ± 0.11) x 10 

3-1-1 1 

cm molecule s , where the error is twice the standard deviation of the 

mean. 

DISCUSSION 

The earlier results for k 1 have been reviewed in the Introduction. The 
present result at 298 K, k 1 = (0.99 ± 0.11) x 10” 11 cm 3 molecule” 1 s” 1 , is in 
excellent agreement with that of Atkinson and Pitts 10 who report (0.94 ± 0.10) 
x 10 1 cm 3 molecule 1 s These workers used the same technique as the 
present (FP-RF). 

Results reported at this temperature from absolute®’^ and relative 1 1 ’ 13 

8 9 11 

DF-MS measurements at low pressure suggest a k^ value > / '30-50/& higher and 

'*' 35 % lower 13 than those obtained by FP-RF. The relative FT1R result of Niki 
et al., 1 ^ which employs an entirely different type of experimental system, is 
in closer agreement with the discharge flow than with the FP-RF results. 

Thus, a disagreement, which may arise from unknown systematic errors, exists 
for this reaction at 298 K. 

The present result, k. = (1.05 ± 0.11) x 10 -11 cm 3 molecule” 1 s” 1 for 228 

10 

< T < 362 K, extends those of Atkinson and Pitts who report (0.94 ± 0.10) x 

10” 11 , (0.94 ± 0.10) x 10" 11 and (1.03 ± 0.11) x 10" 11 cm 3 molecule” 1 s” 1 at 

299, 356, and 426 K, respectively. Although these authors preferred to 

represent their data by the temperature dependent expression, k 1 = 1,25 x 

10” 11 exp(-88 ± 150/T) cm 3 molecule” 1 s” 1 , the individual values differ 

amongst themselves by less than the error at any of their experimental 

temperatures. Thus, a simple average is an adequate representation. If the 

present results and those of Atkinson and Pitts 10 are combined, the rate 

-113 -1 

constant can be best represented by k 1 = (1.01 ± 0.10) x 10 cm molecule 
s 1 for 228 < T < 426 K. This result is to be contrasted to the only other 

temperature dependent study, that of Smith, 13 which leads to k s 6 x 10” 11 

3-1-1 1 

exp(-635 ± 250/T) cm molecule s . This disagreement probably arises, at 

least in part, from the fact that k^ was measured relative to the reaction 

OR + OH + H 2 0 + 0 (4) 

7 
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and the temperatuare dependence of this reaction is virtually unknown over the 

temperataure range of interest here (i.e., 200-500 K) . Although recent 
1 6 

reviews suggest an Arrhenius expression for k ^ consistent with that used by 

Smith, ^ it must be emphasized that these are all based on a measured value 

of k. at 300 K and an estimated activation energy. The older review by Baulch 
** 17 

et al. gives the most detailed summary of the bases for the various 
estimates including their own. 


It therefore appears that reaction ( 1 ) occurs with no energy barrier and 
possibly through two reactive channels, 


OH + H 2 C0 * H 2 0 + HC0 


OH + H 2 C0 ♦ HC00H + H, 


(la) 

(lb) 


reactions (la) and (1b) being exothermic by 3^ and 22 kcal/mole, respectively. 

Path (a) can be considered to be concerted whereas path (b) undoubtedly would 

* 1 6 b 

require the formation of an adduct as in the iaoelectronic case, OH + CgHj., 

forming the radical species, H0CH 2 0. An important question is whether this 

radical is formed with excess vibrational energy (i.e., chemically activated) 

because then the re-dissociation to OH + H 2 C0, in competition with 

dissociation to other products and eollisional stabilization, becomes a 

possibility. Pressure dependence in k 1 might then be possible. The heat of 

* 1 8 

formation of the adduct can be estimated from that of CH 2 (0H) 2 on the 
assumption that the OH bond strength is ^ 105 kcal/mole. Then the calculation 
for the enthalpy change from reactants to adduct yields a value of ^ "22 
kcal/mole. Even though this estimate is uncertain, the adduct, if it is 
formed, could be vibrationally excited. If the energy barrier from adduct to 
I1C00H + H (path(b)) is lower lying than that for re-dissociation to reactants, 
pressure dependence might not be observable. It is also possible that the 
concerted path (a) could arise from such an excited adduct although the direct 
H0 — H — C abstraction path seems more plausible. Clearly the present results 
give no preference as to which reaction channel dominates. 


In an attempt to understand the present experimental results, we have 
carried out theoretical calculations with absolute reaction rate theory where 
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k 1 is given by: 


k 1 = 


kT qt 
h q H 2 CO q C)H 


-E/RT 


(5) 


Our approach in this case is not to use the BEBO method in order to specify 
configuration and energy, since the parameters would have to be adjusted 
(parameterized triplet repulsion) in order to gain the experimental result 
that E = 0. Rather, we ask what do the experimental results imply about the 
activated complex? 


Eq. (5) has been evaluated with known partition functions and the result 
is 


k 1 = (1.70 x 10“ 10 /T 3/2 ) q v t [(ABC)t/(ABC) F ] 1/2 


( 6 ) 


where A, B, and C are the principal moments of inertia for the activated 

t 

complex and formaldehyde, and q y is the vibrational partition function for 
the activated complex (note that q yH CQ and q v , QH are both unity over the 
present temperature range) . Since there is the above noted uncertainty about 
reaction paths, we have considered two models for the activated complex, (1) 

HO — H— CHO where the forming and breaking bonds are col inear, and (2) HO — CH 2 0 
where the forming bond results in a tetrahedral structure. Moments of inertia 
were calculated for these complexes as a Junction of the parameter, r, where r 
is the HO — H bond distance in model (1) and is the HO — C bond distance in 
model (2). The results for model (1) are that [ (ABC)+/(ABC)„] has values 
of 17 > 23, 37, and 56 at r = 1 . 5 » 2.0, 3-0 and 4.0 A, respectively. For the 
tighter model (2), the respective values are 8, 12, 24 , and 40. Thus, for 
either modal with the known value for k., q can be evaluated for the four 

S V 

experimental temperatures of the present study. In most ART calculations two 
low lying bending vibrations in the activated complex make the largest 
contribution to q y . On the assumption that these are degenerate, values can 
be calculated for each model so that agreement with experiment is obtained. 
The calculations are obviously parametric and therefore cannot be used to 
indicate a preference for paths (a) and (b). Briefly, the allowed value for 


the bending vibration is very low for the close .interaction (r = 1.5$, u. is 

— 1 ^ 

55 and 35 cm for models (1) and (2), respectively) and increases to higher 

values for longer interaction (r = 4.0 8. i> b is 115 and 90 cm for models (1) 

and (2),‘ respectively). Since it is generally true that bending force 

constants are greater (and u is larger) for tight complexes, either tight 

complex (r = 1.5 2) would not be indicated. In fact, the exercise suggests 

that in order to tolerate such a large temperature independent rate constant, 

the complex must be quite loose so that "reasonable" values of bending 

frequencies are allowed. 

Lastly, reaction (1) occurs in several systems of environmental 

importance. As mentioned in the Introduction, it plays a significant role in 

high temperature combustion systems and low temperature atmospheric phenomena. 

19 

A recent review of kinetic data for methane-air combustion suggests that k 1 

-11 3-i—i 1 

= 8 x 10 cm molecule s at 1500 K. The increase over the value = 1 
x 10" 11 cm^ molecule”*' s~\ recommended here for 228 < T < 426 K, is probably 
not unreasonable nor unexpected considering the large temperature interval. 

Rate studies over the range 500-1500 K are now feasible using high temperature 
fast flow systems'^ J and would bridge the gap between the presently 
available low and high temperature studies. The impact of the present results 
for low temperature atmospheric modeling is that the rate of tt^CO loss by 
reaction with OH is now shown to be greater than expected from previous 
studies. Thus at 228 K, which corresponds to an altitude of 25 km in 
stratospheric modeling calculations, extrapolation of the Atkinson and Pitts^ 
data leads to an underestimate of k. by about 2056. A far more serious 

13 1 

discrepancy arises from Smith's study since use of the Arrhenius expression 
derived from his data results in a value of k^ at 228 K which is a 'faetor of three 
lower than the present observation at that temperature. 
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FIGURE C APTION 


Figure 1. Examples of reaction-®, and diffusion- H> first order decay plots 
at T = 228 K. The lines are from linear least squares analysis of the 
data points. Both experiments were carried out at a flash energy of 52 


J, a total pressure of 30 torr, and P u n = 60 mtorr. 

h 2 u _ _i 


for the reaction experiment. 

-1 


k . , = 487 + 26 s 

observed 


p h 2 co = °- 94 mtorr 

and k d = 30 + 2 
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Table 1. 

Rate data for 

the flash photolysis-resonance fluorescence study of the 
reaction OH + H 2 C0 

c* 

T 

[H CO] 

[HO] 

[Ar] 

Flash energy a 

No. of 

k 1 b 

- 

(K) 

(mtorr) 

(mtorr) 

(torr) 

(J) 

expts. 

(10“ 12 cnr 5 molecule 


228 

0.50 

40 

20 

24 - 35 

5 

11.23 + 1.56 



0.63 

40 

20 

47 - 52 

3 

11.30 + 0.17 



0.94 

60 

30 

52 

1 

11.60 

! 


1.00 

80 

40 

11 - 35 

5 

11.09 + 0.?7 
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11.22 + 0.98^ 


257 

0.63 

40 

20 

43 - 57 

3 

10.20 + 1.16 

1 


0.76 

40 

20 

47,52 

2 

10.13 + 1.46 



0 . 94 

60 

30 

39,47 

2 

10.17 + 0.22 

> 


1. 14 

60 

30 

43 

2 

9.48 + 0.48 



1.25 

80 

40 

43 

1 

10.48 



1.52 

80 

40 

43 

2 

11.38 + 0.26 
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10.28 + 0. 90 u 


298 

1.00 

40 

20 

43,59 

2 

9.94 + 0.74 



1.13 

60 

30 

43,55 

2 

9.42 + 0.48 



1.50 

60 

30 

36 - 55 

5 

10.28 + 0.93 

* 


1.00 

80 

40 

43,57 

2 

8.12 + 0.70 

} 


2.00 

80 

40 

21 - 57 

3 

10.08 + 0.94 

f 


1.69 

90 

45 

43 

1 

8.74 

s 


1.50 

120 

60 

21 - 58 

4 

9.23 + 1.08 

r 


2.25 

120 

60 

31 - 52 

4 

10.28 + 1.34 



2.00 

160 

80 

21 - 57 

4 

10.93 + 0.72 

k 






25 

9.86 + 1.15° 


362 

0.62 

20 

20 

24,39 

2 

8.55 + 0.56 



0.94 

40 

20 

43,52 

2 

11.29 + 1.22 

i 


1.23 

40 

20 

24,39 

2 

10.51 + 3-25 



0-92 

30 

30 

24,39 

2 

11.58 + 1.44 

1 


1.41 

60 

30 

43,52 

2 

10.09 + 2.28 

. 


1.85 

60 

30 

24,39 

2 

10.66 + 0.89 

r 

1 


1.23 

40 

40 

24 - 43 

4 

10.50 + 0.38 

a 


1.88 

80 

40 

43,52 

2 

9.58 l 2.03 

t 


2.47 

80 

40 

24,39 

2 

11-33 + 1.81 



i , a 

' b 157.5 nm filter on flash lamp 

§•;. c Error limit in is one standard deviation 

F- Mean value of k at that temperature 

* * 
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